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■ Abstract 

pL ; 

We study an asymptotic behaviour of the principal eigenvalue for an elliptic operator 
. with large advection which is given by a gradient of a potential function. It is shown that 

the principal eigenvalue decays exponentially under the velocity potential well condition 
as the parameter tends to infinity. We reveal that the depth of the potential well plays 
Oh ' an important role in the estimate. Particularly, in one dimensional case, we give a much 

. more elaborate characterization for the eigenvalue. Some numerical examples obtained by a 

^ ' characteristic-curve finite element method are also shown. 

» 

<^ ■ 

1 Introduction 

The following elliptic eigenvalue problem with a large drift term is considered in this paper. Let 
CN ■ Q be a bounded Lipschitz domain in (n > 1). For a given vector field a = (oi, • • • , an)"^ S 

^ L°° {fl,W^), we consider an elliptic eigenvalue problem with a parameter p E M : 

-Au{x) -\-pa{x) ■ Vu{x) = Xu{x) {x € 0,) 

(1.1) 

g ; [ u{x) = {x€ dn) 

00 

Q I We assume that the eigenfunction u {u ^ 0) belongs to the Sobolev space Hq{Q,). Then, from 

the elliptic regularity theorem (see [9] for example), u S W^^^{U) holds for arbitrary s E [l,oo). 
- Although A and u{x) are complex-valued in general, we can define a positive real-valued 

r> I principal eigenvalue A = \i{p) and the corresponding real- valued eigenfunction u = ui{-,p), 

' which is uniquely determined by the condition: 

u{x) > {x £ il), maxu(a;) = 1 . (1.2) 

Some related results for the principal eigenvalue and eigenfunction will be collected in Theo- 
rem 13.11 and Theorem 13.71 A detailed and systematic study for the principal eigenvalue for 
general second order elliptic operators in general domains can be found in Berestycki, Nirenberg 
and Varadhan [3|. 

As we will see in Section [2] through several numerical examples, the principal eigenvalue 
is closely related to the decay rate of the solution of the corresponding nonstationary linear 
advection-diffusion equation (12. Ij) . 

Our interest in this study is asymptotic behaviours of \i{p) as p ^ oo. Such large effects 
from advection or drift term under small diffusion appear in many actual physical problems and 
often cause some difficulties in their numerical simulations and analysis. 
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This phenomenon was first studied by Ventcel' and Friedman [8j with Ventsel'-FreidUn's 
probabihstic approach [19]. In [8], the fohowing result is obtained. If O is a bounded domain of 
C^-class, and if a belongs to C^{Q) and satisfies the condition: 

a{x)-v{x)>0 {xGdi}), (1.3) 

where v denotes the outward unit normal vector on dQ, there exist ci and C2 such that 

< ci < lim inf - log ^ , , < lim sup - log , , . < C2 . 
p^oo p \i[p) p^oo P MkP) 

This estimate is equivalent to 

> 0, ^po e M s.t. e-(^2+e)p < x^{p) < e-(=i-^)P > Po), 

which means that the principal eigenvalue becomes exponentially small under the condition 

(USD. 

In one dimensional case, asymptotic behaviours of kth. eigenvalues Xk{p) {h G N) are studied 
in |5| by means of matched asymptotic expansions. In connection with a singular limit analysis 
of Sturm-Liouville two points boundary value problems, precise approximations of eigenvalues 
and eigenfunctions including exponentially small principal eigenvalue are obtained there. 

For a given vector field a G L°°(r2,M"), if the principal eigenvalue \i{p) satisfies the condi- 
tion: 

and ^po G ^ s.t. < \i{p) < e^^^'P i^p > po), (1.4) 
or its equivalent condition: 

lim inf - log -r—r^ > 0, 

p-^oo p M[p) 

then we call it exponential decay phenomenon of the principal eigenvalue. In addition to some 
numerical examples in Section [2l we will give a biological interpretation in a chemotaxis model 
for more physical pictures of the exponential decay phenomenon. 

Besides the exponential decay phenomenon, the principal eigenvalue Xi{p) exhibits various 
asymptotic behaviour as p — > oo. Devinatz, Ellis and Friedman [5] investigated inequalities 
type arguments and maximum principle type arguments, and they obtained some estimates of 
\i{p) from above or from below by a term Cp"' with 7 G (0, 2]. See also jl8j for related results. 

In [2j, Berestycki, Hamel and Nadirashvili proved that, if dQ is of C^-class and a G L°°{Q, R") 
satisfies diva = in P'($7), then limsupp_,oo Ai(p) < 00 if and only if there exists ijj G H^{^1) 
such that V' ^ const and a ■ ViJj = a.e. in 0,. 

In our study, we focus on the exponential decay phenomenon and give some new estimates 
for it. We assume existence of a velocity potential of a (condition (j3.1|) ). which allows us to 
use the Rayleigh quotient and inequalities type arguments. One of our main results shows 
us that existence of a potential well (Definition 14. 6p implies the exponential decay phenomenon 
of the principal eigenvalue, and that the depth of the potential well gives an estimate of the 
constant c in ()1.4p from below (Theorem 14. 7p . This observation is extended to the exponential 
decay phenomenon of the mth eigenvalue (Theorem 14. Sp and to a precise asymptotic behaviour 
of Xi{p) in one dimensional case (Theorem 15. 2p . 

The organization of this paper is as follows. In Section [2 we show some numerical examples 
of the exponential decay phenomenon in two dimensional case. These numerical profiles of the 
principal eigenfunctions ui{x,p) will be helpful in our analysis later. We collect some funda- 
mental facts on the principal eigenvalue in Section [3l In Section HI we prove several estimates 
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for asymptotic behaviours of Ai(p) including an alternative simpler proof of the exponential 
decay phenomenon. In Section O we give more precise asymptotic behaviour of Xi{p) in one 
dimensional case, under a different assumption than one of [3]. Additionally as an application, 
we give a biological interpretation of the exponential decay phenomena, life span of a biological 
colony under the chemotaxis effect in the last section. 

2 Numerical examples of exponential decay phenomena 

For the principal eigenvalues Xi{p) which is defined in (jl.ip . we give several examples of typical 
asymptotic behaviours of Ai(p) as p ^ oo for fixed vector fields a{x), particularly focusing our 
interests on the exponential decay phenomenon. We start from the simplest case that a is a 
constant vector field. 

Example 2.1. Let a £ M" be a constant vector. Then we can easily check that 

lap p 
Xi{p) = Ai(0) +/^, ui{x,p) = e2''-''ui{x,0), 

where Ai(0) and ui{x,0) are the principal eigenvalue of —A with the Dirichlet boundary con- 
dition on Q and its corresponding eigenfunction. Moreover, not only the principal eigenvalue 
but also all eigenvalues are exactly shifted by (|ap/4)p^, i.e., Xk{p) = Afe(O) + (|ap/4)p^ for all 
k £ 'N. In this case, the asymptotic behaviour of Xi{p) (or Ajt(p)) is of 0{p^) as p ^ oo. We will 
see in the next section that this type of 0{p^) behaviour is most common, for example if the 
vector field a{x) has no singular point. Of course, there is no chance to have the exponential 
decay phenomenon of principal eigenvalues as p — > oo. 

On the other hand, it seems to be difficult to give an example of the exponential decay 
phenomenon similarly only by using elementary function. 

In this Section, we give some examples of the exponential decay phenomenon with the help 
of numerical simulations. For this purpose, we compute the following initial-boundary value 
problem of parabolic type: 

ut- l^u + p a{x) -Vu = (a; G il, < t < oo) 
< u{x,t) = {xedn,t>0) (2.1) 

u{x,0)=uq{x) (a; G ri), 

where uq is a given positive initial function. In the following numerical examples, the domain 
is chosen as Q = (—1, 1) x (—1, 1). It is well-known that u{x,t) uniformly converges to zero as 
t — > oo and its asymptotic behaviour is more precisely given as 

u{x,t) Ce-^'^P^^ui{x,p) ast^oo, (2.2) 

where the constant C is positive if U(){x) > for x £ Q. 

Ifp >> 1, since (j2.ip is related to exponentially small positive eigenvalue and is a convection- 
dominant problem, it requires us to adopt some upwinding technique for its reliable and reason- 
able numerical simulation. The following simulations were computed by means of a characteristic- 
curve finite element scheme (103) in Chap. 3 of [14] with piecewise linear triangular elements. 
This is a upwind implicit scheme of first order based on the characteristic curve approximation: 

/ . \ , . „ , , s u(x, kr) — u(x — pa(x)T, (k — 1)t) ^, , 
ut{x, kr) +pa{x) ■ Vu{x, kr) = ^ ^ \ i ^\ _^ q(^^^^ 

T 
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where r > is a small time increment. For more details and the stability analysis, see |13] and 
[H] etc. 

We use the following vector field with compact support on M^; 



We remark that a{-;R) £ VF1'°°(M2; M^) and ||q;(-, i?)||ioo(iR2) = 1. 

Example 2.2. We set a{x) := a{x; 1/2) and uq = 1. For p = 0, 10, 20, 30, 40, we numerically 
computed (j2.ip in the time interval < t < 1.0 and the results are shown in Figures [THlll For 
each p, time evolution of u{x,t) and a profile of principal eigenfunction ui{x,p) are shown in 
the figures. 

These finite element simulations were computed with time increment r = 1/2000 on a 
unstructured triangular mesh of $7 which consists of 11292 total nodal points (degree of freedom 
including boundary points) and 22182 triangular elements with the average edge length 0.021. 
For the mesh generation, we used FreeFem++ [7]. 

Time evolutions of u{x, t) are drawn with fixed X2 = 0, where the horizontal axis stands for 
— 1 < xi < 1. In each simulation, for t > 0.2, u{x,t) decreases monotonically and keeps almost 
same profile, i.e., u{x,t) exhibits the asymptotical form (j2.2p . Comparing these figures, we can 
see that the decay speed of u{x, t) becomes extremely slow in cases oi p = 20, 30, 40. Especially 
in case of p = 40 (Figure [9|), it seems to be almost stationary after t = 0.2 in the figure. This is 
a typical example of the exponential decay phenomenon. 

The principal eigenfunctions drawn in the figures are normalized profiles of u{x, t) at t = 1.0, 
which are considered as profiles of principal eigenfunctions ui{x,p). The sections of normalized 
ui{x,p) on the line X2 = are drawn in Figure [TTJ We remark that, in case of p = 0, the 
principal eigenvalue and eigenfunction are given by 



The support of the flow field consists of two disjoint sets {|a; — < 2/5} and — a;2| < 1/4}. 

For p = 0, 10, 20, • • • , 100, we numerically computed ()2.1I) in the time interval < t < 1.0 
and the results are shown in Figures [T214221 For each p = 10, 20, 30, 50, 100, time evolution 
of u{x,t) and a profile of principal eigenfunction ui{x,p) are shown in the figures. These 
finite element simulations were computed with same time increment and triangulation of 17 as 
Example 

Time evolutions of u{x,t) are drawn on the line 6xi — 10x2 + 1 = through xi and X2, 
where the horizontal axis stands for — 1 < a^i < 1. Similarly to the previous example, this also 
exhibits the exponential decay phenomenon. 

The principal eigenfunctions drawn in the figures are normalized profiles of u(x, t) at t = 1.0, 
which are considered as profiles of principal eigenfunctions ui{x,p). The sections of normalized 
ui{x,p) on the line 6xi — 10x2 + 1 = are drawn in Figure [22l 

From this and previous examples, the profiles of the eigenfunction ui{x,p) seems to converge 
to a limit profile as p ^ oo. In particular, the limit profile seems to be fiat on the support of a 
in both cases. 




(0 < \x\ < R) 



{x = 0, or \x\ > R) 



^i(O) = y, 'ui(a;,0) = cos — cos— . 
Example 2.3. We set Xi := (1/2, 2/5) and X2 := (-2/3, -3/10) and define 

a{x) := a.{x — Xi; 2/5) + 2a{x — X2; 1/4). 
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t=0.0 




Figure 1: Time evolution of Example 12.21 Figure 2: Eigenfunction of Example 12.21 for 

on the line X2 = for ]5 = 0. p = 0. 



t=0.0 




Figure 3: Time evolution of Example 12.21 Figure 4: Eigenfunction of Example 12.21 for 

on the line X2 = for p = 10. p = 10. 
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Figure 7: Time evolution of Example 12.21 
on the line X2 = for p = 30. 




Figure 8: Eigenfunction of Example 12.21 for 
p = 30. 




Figure 9: Time evolution of Example 12.21 Figure 10: Eigenfunction of Example 12.21 

on the line X2 = for p = 40. for p = 40. 
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Figure 11: Eigenfunctions of Example 1 2 . 2 1 on the line 
X2 = 0foTp = 0, 10, 20, 30, 40. 
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Figure 12: Time evolution of Example 
on a line for p = 10. 




Figure 13: Eigenfunction of Example 
for p = 10. 




Figure 14: Time evolution of Example 
on a line for p = 20. 




Figure 15: Eigenfunction of Example 
for p = 20. 




Figure 16: Time evolution of Example 
on a line for p = 30. 




Figure 17: Eigenfunction of Example 
for p = 30. 
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t=0.0 




Figure 18: Time evolution of Example 12.31 Figure 19: Eigenfunction of Example 12.31 

on a line for p = 50. for p = 50. 




Figure 22: Eigenfunctions of Example 12.31 on a line 
for p = 0, 10, 20, • • • , 100. 



3 Fundamental tools for eigenvalue problems 

In this section, we collect several fundamental facts and tools for the eigenvalue problem. For 
simple notation, we fix p = 1 without loss of generality throughout this section. 

Theorem 3.1. In a bounded Lipschitz domain Q C (n > 1), we assume a G L°°(r2,M"). 
Then, there uniquely exists a positive number Ai and ui £ H^^^{Q) n {{^{fl) such that A = Ai 
and u = ui satisfy (II. ip with p = 1 and (II. 2p . Moreover, Ai is given by the min-max formula: 

Aip{x) + a{x) ■ Vip{x) 

Ai = max mf ' 



where max^>o is taken over all positive ip S W^^^{Q). 

We remark that inf and sup always denote essential infimum and supremum in this paper. 
For the proof of the theorem, see Section 2.8 of |15j and [3]. Since Ai = min^ ReA holds for 
any complex eigenvalue A of —A + a • V, Ai and ui are called the principal eigenvalue and 
eigenf unction. The next corollary immediately follows from Theorem 13.11 

Corollary 3.2. Under the condition of Theorem \3.1l we consider a Lipschitz subdomain 0,' C 0, 
and define X'l > as the principal eigenvalue for — Au+a-Vu = X'lU in 0,' with the zero Dirichlet 
boundary condition on dO,' . Then Ai < A'^^ holds. 

Throughout the following sections, we assume the existence of a velocity potential of a: 

^beW^'^^in) s.t. a{x) = Vb{x) (a; e 17). (3.1) 

Since the existence of a velocity potential is not essential for the exponential decay phenomena 
which is our aim in this paper (see [8j, for example), the assumption ()3.1|) is rather technical 
but allows us energy-based arguments with the help of the following Liouville transform. 

Under the condition ()3.ip and for q € L°°(J7), we consider the following three elliptic eigen- 
value problems with the Dirichlet boundary condition. 



—Au{x) + a{x) ■ Vu{x) = Xu{x) {x G Vl) 

u{x) = {xe dn) 

—Av{x) — dw{v{x)a{x)) = Xv{x) {x G Q) 

v{x) = (xe on) 

—Aw{x) + q{x)w{x) = Xw{x) {x e 17) 

w{x) = {xe on) 



(3.2) 



(3.3) 



(3.4) 



Proposition 3.3 (Liouville transform). Suppose the condition (|3.ip . Then the above three 
eigenvalue problems are equivalent to each other under the following relations: 

b b 11 

e2v = e~2u = w, q = —-Ab + -\Vb\'^ . 

We omit the proof, since it is shown by straightforward substitution. We remark that the same 
Liouville transform is valid even for the time dependent problem ()2.ip . 
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For the self-adjoint eigenvalue problem ()3.4p with q G L°°(r2), there is a well-known min-max 
principle for the Rayleigh quotient. We define 

:= -A + q, Dom(L,) = HI{Q), (3.5) 
H-^[n){L,w,w)Hl(n) _ !^{\Vw{x)\'' + q{x)w{xf)dx . _ ^i.^.. 

Jq[w) .= - J. [w£HQ[il)). 

{w,w)L2(n) J^w{xydx 

Theorem 3.4. Forq G L'^{Q), there exists a complete orthonormal system of L'^{Q), {wk}k<^n C 
Hl{VL) n Hf^^{VL) and {\k}km C M (Ai < A2 < • • • ) such that 

LqWk = XkWk (A; G N). 

The kth eigenvalue Xk is characterized by the min-max formula: 

Xk = min max Ja(w), 

dimX=kweX\{0} ^ 

where the minimum is taken over whole k-dimensional subspace X in Hq{Q,). In particular, the 
principal eigenvalue Ai is simple and given by 

Ai = min Jq{w). 

«;6Hi(n)\{0} 

Moreover, an eigenfunction w G Hq{^) n ffj^^(il) corresponds to the principal eigenvalue if and 
only if w has no sign change in Vt. 

For the proof of this theorem, see Section 8.12 of [9j and pLj etc. 

From the last assertion of this theorem, without loss of generality, we can assume that 

wi{x) > (a; G 0). 

Theorem 3.5 (comparison theorem for principal eigenvalues). For q, q G L°°{fl), the principal 
eigenvalues of Lg and Lg are denoted by Ai and Xi, respectively. Then 

mi{q - q) < Xi-Xi < sup(g - q). 
" n 

In particular, if q{x) > q{x) {x G il) and \\q — q\\L^(^[i) 7^ 0, then Ai > Ai. 

Proof. Let wi and wi denote the eigenfunctions corresponding to Ai and Ai, respectively. We 
have 

Al = min Jq{w) < Jq{wi) = Jq{wi) + {Jq{wi) - Jq{wi)) 

T ^ IniQjx) - q{x))wi{x)^dx ~ 
= + r wUxWx <Ai+sup((z-g). 

In the same way, we also have 

Al < Al -I- sup{q - q) = Xi- inf(g - q). (3.6) 

n ^ 

These inequalities imply the first assertion. Furthermore, q > q and q ^ q then, from the 
positivity of wi and wi, we have 

Xi{wi,Wi)l2^q) = {LqWl,Wl)L2(^^) = {wi , LgWl) l2(^^-) > {wi,LqWl)L2(^^) = Al (uil , 7Z;i)i2(Q) , 
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and Ai > Ai follows. □ 

We define the elliptic operator Kb by 

i^fe := -A + V6 • V, Bom{Kb) = H^{n), 

and denote by its adjoint operator with respect to L'^{Q): 

K; := -A - div(-V6), Dom{Kl) = H^{n). 

Although the following results concerning the operator will not be used in our analysis, we 
include them below for our systematic description. Another reason to include them is that the 
eigenvalue problem ()3.3p corresponding to is related to various important applications as 
well as (13. 2p for Kb- An example of such applications will be shown in Section [H 
We introduce the following weighted inner product of L^(r2); 

(f,n)i2(n) := / e^^'^^v{x)u{x)dx, {v, ueL^{Q)), 



and we define Lf^{^) '■= (L^(r2), (•, Ol^^q)) which denotes a Hilbert space L?'{Q.) with this inner 
product. We also define 

Le^^'^^UVvix^^^ - J\bix)vixf}dx 
Then we have the following proposition. 

Proposition 3.6. We suppose that u, v, ip £ Hq{Q) n H'^{Q) (u, v ^ 0). Then we have 
{KbU,ip)L2 = {u,Kb'p)L2 tn), Ib[u) = — ^ , 



Namely, Kb and K^ are self adjoint in L^^(r2) and respectively, and their Rayleigh quo- 

tients are given by Ib{u) and II {v). 

Proof. For ip = u or v,hy the integration by parts, we obtain 

/ e'^^\V^\^dx = - [ div(e'^''V^)^dx= [ e'^'' {-Aip ±Vb ■ V^) ipdx. 
Hence, we have 

J^e-\-Au + Vb-Vu)udx _ J^e-^iKbu) udx _ 
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_ e^K-A^; -Vb-Vv)v- (Afe) v^}dx _ J^e\Klv)vdx _ iK;v,v)q(^^) 

The self-adjointness of Kf, in L^^(r2) and the one of in L^(r2) are also checked by direct 
calculations. □ 

We have the following characterizations of the eigenvalues of ()3.2p and ()3.3p . 

Theorem 3.7. For b G M^^'°°(r2), ^/lere exist {AfejfcgN C M (Ai < A2 < • • • ); '^'^^ {^^fcjfceN; 
{^JfcjfceN C -H'o(^) n ^^ioc(^) s^^c/i that 

KbUk = XkUk, Klvk = \kVk (A;GN), 

and that {uk}k&] and {vk}k<^N ctre complete orthonormal systems of LPi^{rL) and Ll{Q.), respec- 
tively. The kth common eigenvalue Xk is characterized by the min-max formula: 

Afc = min max Ih(u) = min max Ih(v), 

dimX=kueX\{0} dimX=fct,6X\{0} 

where the minimum is taken over whole k-dimensional subspace X in Hq(Q). In particular, the 
common principal eigenvalue Ai is simple and given by 

Ai = min hiu) = min Ibi^)- 
«eHi{Q)\{o} veH^{n)\{o} 

Moreover, an eigenfunction u for Kf, {v for K^) corresponds to the principal eigenvalue if and 
only if u (v) has no sign change in ft. 

Proof. The assertions follow from Proposition 13.31 and Theorem 13.41 □ 

Prom the last assertion of this theorem, without loss of generality, we can assume that 

ui{x) > 0, vi{x) > (xefl). 

As a last remark in this section, we introduce the following proposition. We also omit the proof 
since it is shown by a direct calculation. 

Proposition 3.8. Let Q := H^Lj^Tj C M", where li is a bounded open interval for i = 1, • • • ,n, 
and let a{x) = (ai(xi), 02(2:2), •• • ,an{xn))'^ for x = (xi,-- - ,a;„)'^ G fl with ai G L'^{Ii). 
Then the principal eigenvalue Ai and the eigenfunction ui{x) of (j3.2p are given by the following 
formula: 

n 

U ' 

i=l 1=1 



Ai=^A^*\ ui{x) := Yl^iixi) (xeQ), 



where a[*^ and ipi are defined by 



(x) + ai{x)ip'i{x) = X'i^ipi{x) {x G Ij) 
ipi{x) =0 (x G dli) 

, ipi{x) > (x G li). 
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4 Asymptotic behaviour of principal eigenvalues 

In this section, we consider the singular perturbation problem of the principal eigenvalues (11.11) 
under the velocity potential condition (j3.ip . Henceforth we assume that p is a positive parameter. 

Due to Proposition 13.31 (|l-ip is equivalent to each of the following three eigenvalue problems 
under the zero Dirichlet boundary condition: 

KpbU = — An + pVb ■ Vu = Xu in fi, 



K*i^v = —Av — pdiv{vVb) = \v in Q, 
Lq(p)W = — Att) + q{p)w = Xw in Q, 

where we define 

2 

q{x,p) := -|diva(cc) + ^|a(a;)|2, (4.1) 

and we abbreviate q{-,p){G L°°{Q)) as q{p). 

Under the condition (j3.ip . the kth eigenvalue (their multiplicities are counted) of these 
equivalent eigenvalue problems is denoted by Xkip) for a parameter p > 0. 

We start from the following simple consequence of the comparison theorem. 

Theorem 4.1. Let Xq > be the principal eigenvalue o/— A/?, which is Laplacian with the zero 
Dirichlet boundary condition. Under the condition (13. ID . following two estimates hold: 

p 

Ai(p) > - - sup(diva(a3)) {p > 0). 

1 ■ r I / m2 ^ r ■ f Al(p) , Al(p) 1 I . m2 

- mi a ic) < limmf — tt— < limsup — tt- < - sup auc) . 

4 teen' p^oo p2 p2 4 3,gn 

Proof. Applying Theorem 13.51 to -^^g(p) and Lq = — A/j, we have 

inf g(a;,p) < Ai(p) - Aq < supg(a;,p), (4.2) 

and 

/ 2 \ 2 

Ai(p) > An + inf ( -|diva + ^\a\'^ J > Ao - | sup(diva) + ^ inf jap, (4.3) 

--diva + y |ap j < An - - inf(diva) + — sup |ap. (4.4) 

The first assertion of the theorem follows from (j4.3p . Dividing (|4.3p and (j4.4p by and taking 
the limit p ^ oo, we can also derive the second assertion. □ 

Corollary 4.2. There is no exponential decay phenomenon ifmf\a{x)\ > or supdiva(a;) < 0. 

The estimate (14. 2p can be improved as follows. 
Lemma 4.3. For pi > P2 ^ 0, we have 

Pi + P2 

inf q{x,pi +P2) < (Ai(pi) - Ai(p2)) < sup q{x,pi +P2). (4.5) 

xen P1-P2 xen 
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Proof. From the definition of q{x,p) ()4.ip . the equahty 

q{x,px) - q{x,P2) = : q{x,pi +P2), 

Pi +P2 

holds. Hence, the assertion follows from Theorem 13.51 □ 
Prom (j4.2p and (|4.5p . it is natural to consider the following condition for a: 
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'po > s.t. inf^q{x,po) > 0. (4.6) 
Since 

^ = --diva + -|a| , 

the inequality 

^±2Pll>^±lP^ (m>P2>0), (4.7) 

Pi P2 

holds. Hence, the condition (|4.6p implies mixenQix,p) > for all p > pQ. From Lemma 14.31 
we have the following theorem. 

Theorem 4.4. Suppose the condition (j4.6p . Then, \i{p) is nondecreasing for p > pq/2. In 
particular, the exponential decay phenomenon does not occur in this case. 

Proof. We assume that pi > P2 ^ Po/2. Then we have pi + P2 > Po- From the first inequality 
of Lemma 14^ and (14.71). we obtain 



m{Pi) - m{P2) > : mf g(pi + P2) = [Pi - P2) mf ■ > [pi -p2)int > 0. 

Pi +P2 Pi+ P2 Po 

□ 

We notice the following necessary condition for (j4.6p . 
Proposition 4.5. The condition (|4.6p implies that 

min b(x) = min b(x) (0,': an arbitrary subdomain of Q). (4.8) 

xeW xedQ' 



Proof. We remark that the condition ()4.6p is equivalent to JCb > 0, where /C is a linear differential 
operator defined by 

/Cm := -An + yV6- Vn. 

Since 6 is a supersolution with respect to fC, from a consequence of the weak maximum principle 
([9], Theorem 3.1), the condition ([^ follows. □ 



An inverse condition of ()4.8p is given by the following potential well condition. 

Definition 4.6. For a fixed velocity potential b £ W'^'°°{Q), a Lipschitz (nonempty) subdomain 
ri' C O is called a potential well if the condition 

min b(x) < min b(x), 
xeW xedn' 



is satisfied. Furthermore, 
is called the depth of a potential well W 



bo := min b(x) — min b(x) > 0. 
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We remark that if there exists a potential well, the condition ()4.8p does not hold, and neither 
does (|4.6p . Moreover, according to the next theorem, existence of a potential well implies the 
exponential decay phenomenon of principal eigenvalues. 

Let Q' be a subdomain of C^-class and let u denote the outward unit normal vector on 
If it satisfies the condition: 

a{x) ■ iy{x) > {xe dQ'), (4.9) 

which is a sufficient condition for the exponential decay phenomenon in W obtained in [8], then 
the exponential decay phenomenon occurs also in 0, due to Corollary 13.21 Actually, under the 
condition (j3.ip . the condition (j4.9p implies that J7' is a potential well. 

Theorem 4.7. We suppose that there exists a potential well il.' with depth bo > 0. Then an 
exponential decay phenomenon occurs: 

lim inf - log — rr ^ ^0 • 

p^oo p Xi{p) 

In other words, for any oo £ (0, 6o), there exists C > such that 

< Xi{p) < Ce-'^P (p>0). 

Proof. Without loss of generality, we can assume mm^^jyb(x) = 0. From Theorem 13.71 Xi{p) 
is given by 

Xi{p) = min Ipb{u). 
Hence, for arbitrary u € Hq{^) {u ^ 0) and /? G M, we have 

Let us choose (3 satisfying 0</3</3 + a;<6o- Since b G C^{Q) and ming^' b = bQ > (3 + lj, 
there exists e > such that b{x) > /? + for all a; G N'^{dQ,'), where 

N%dn') := {x G n'; dist(a;, < e} , dist(a;, := min{|a; y G dQ'}. 

We define u G H^{n) n C^^^{Tl) by 

'0 {xen\n') 

u{x):={ e-^dist{x,dn') [x € N'{dQ')) (4.11) 

{x G n'\N'{dn')). 



and substitute it to (fiJOj) . Since \Vu\ = a.e. in N^{dn') and \Vu\ = in S7 \ N'^idn'), we 
obtain 

" In'\NHdn') e-P^^(-)-^)dx " \{x G f)'; 6(a.) < f)}\ 

where 

C := £ f]'. < |iv^(ar?')| . 

□ 

The above theorem can be extended to first m eigenvalues Ai(p), • • • , Am(^'), if there are m 
potential wells. 
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Theorem 4.8. Let C (j = 1, • • • , m) be disjoint potential wells (QinO,j = for i ^ j), and 
let loq be the minimum depth of the potential wells ^lj {j = 1, - ■ ■ ,m). Then, for any uj G (0, ujo), 
there exists C > such that 

< Ai(p) < A2(p) < • • • < \m{p) < Ce~^P {p > 0). 

Proof. For each j = 1, - ■ ■ , m, we define Uj S Hq{Q) n C^'^{Q) with supp{uj) = Qj in the same 
manner as (j4.1ip . Let X be a linear subspace of Hq{Q) generated by functions ui, - ■ ■ , Um- It is 
easy to see dimX = m. Due to the min-max formula in Theorem 13 .71 the mth eigenvalue \m{p) 
is estimated from above as 

Am(p) < max Ipb{u) = max Ipb{ciUi ^ h CmUm), (4.12) 

«ex\{o} cj 

where in maxc^. , the maximum is taken over all (ci, • • • ,Cm,) G M™ \ {(0, • • • ,0)}. Since Vtj are 
disjoint, we have 



Iphiciui -\ h 



( I yc'\u,{x)A dx E^i / e-^'^^¥,i^)\'dx 

Jn / j=i •'^^ 



e-pb{^)|Vu^.(a;)|2(ia; 



3=1 ^ I e-P^^'^\uj{x)\'^dx I e-P^^'^'^\yuj{x)\'^dx 



< max 



Yd] / e~^'^-\u,{x)\^dx 



where we have defined dj := Cj{J^ e p^^'^^ \uj {x)\'^ dx)^/"^ . 

Similar to the proof of Theorem 14.71 choosing sufficiently small e > 0, we can show that the 
last term is bounded by Ce~'^P from above. □ 



5 Precise asymptotic behaviour in one dimensional case 

In this section, we introduce a more precise estimate for the exponential decay phenomenon of 
the principal eigenvalues in one dimensional case. In [12], some of results in this section have 
been already announced in Japanese by one of the authors. 

We consider the following one dimensional eigenvalue problem on Q = (—1,1), where I > 0. 

—u"{x) +pa{x)u'{x) = \u{x) {—I < X < I) 
< u{-l)=u{l) = (5.1) 

u{x) > {-I < X < I). 

The principal eigenvalue is denoted by Xi{p) > 0. In the following argument, we fix one of the 
primitive functions of a G L°°{-l,l) and denote it by 6 G VF^'°°(-Z,/) = C°'^([-Z, /]), i.e., 

b'{x) = a{x) a.e. x £ (—1,1). 
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We also define 



bo '■= ^2 — bi, bi := min b(x), 62 := max b(x), 

0<x<l 0<x<l 



Bi:={x£ [0,1]; b{x) = h} (i = 1, 2). 
We assume the following assumption. 

Assumption 5.1. The function a{x) is an odd function belonging to L°°{—l,l), and it satisfies 

max X < min y . 

xeBi J/6-B2 

A typical example of a(x) satisfying Assumption 15 . 1 1 is 

a{x) = ^ b{x) = -1x1", (5.2) 
X a 

for a > 1 and I > 0. In case that 

a{x) = sinx, b{x) = —cosx, (5-3) 
it satisfies Assumption 1 5 . 1 1 if < / < 27r. Another example is 

a{x) = x^-x, b{x) = ^{x^ - if , (5.4) 

which satisfies Assumption 15.11 if / > \/2- 

In the following arguments, A{p) ~ B(p) as p — > cxd means that A{p) / B[p) ^ 1 as p ^ cxd. 
Now we present the main result of this section. 

Theorem 5.2. Under Assumption \5.1l we have 

Ai(p) ~ (^^'e"P^(^)(ix^ (^^'eP^(^)dy^ as p ^ 00. (5.5) 

Proof. We first remark that the principal eigenfunction ui{x) defined by (j5.ip is an even function 
since a{x) is odd. By means of the change of variables from x G [0, to r e [0, 1] by r = p{x): 



p{x):=- r eP^'^yUy {0<x<l), a := [ e^^^^^dx, 

0- Jo Jo 

Ui{r) := u^{p-\r)) (0 < r < 1), 

it follows that Ui{r) and \i{p) satisfy 

' - U'{{r) = \i{p) z{r) Ui{r) (0 < r < 1) 

C/{(0) = 0, [/i(l)=0 (5.6) 

t Ui{r) > (0 < r < 1), 
where 

z{r) := a^e-^PK/'-'M) > (0 < r < 1). 
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We remark that the quantity in the right hand side of ()5.5I) is given by the following equality 



I z{r)dr = (^j^ e-P^^^^dx^ (^j^ eP^^y^y 



which can be checked by direct calculation. 

The eigenvalue problem (j5.6p is selfadjoint and it is found that \i{p) is given by minimum 
of the Rayleigh quotient: 

, , , ^.Wm'dr . !',y{rrdr 

Xi{p) = = mm , (5.7) 

/q z{r)Ui{rYdr v'ey z(r)ip{r)'^dr 

where :={(/? G i?^(0,l); 99(1) = 0, ip^O}. 
Since 

U[{r) = U[{0)+ [/('(s)ds = -Ai(p) / z{s)Uiis)ds < (0 < r < 1), 
Jo Jo 

we have maxC/i = Ui{0). Without loss of generality, we can assume that maxUi = Ui{0) = 1. 
Among functions if S H^{0,1) with (p{0) = 1 and (p{l) = 0, the minimum of the Dirichlet 
integral \ip'{r)\'^dr is achieved by (p{r) = 1 — r. Hence, from the first equality of (j5.7|) . we get 
an estimate from below: 

/q z{r)Ui{rYdr z{r)dr z{r)dr 

On the other hand, from the second equality of (j5.7p . substituting ip{r) = 1 — r, we obtain 
an estimate from above: 

. jjWirrdr ^ /J 1(1 -r)fdr 1 
Xi{p) = mm^^^ < — p = — . (5.9) 

vsV z{r)ip{r)'^ dr •2(r)(l — r)'^dr z(r)(l — r)'^dr 

From these two estimates ()5.8p and ()5.9p , if we can show that 

(} z{r){l — r^dr 

^° \, \ ^1 as p^oo, (5.10) 

Jo ^y'^jdr 

then we obtain the assertion of the theorem. 
Let us show (j5.10p . We define 

Z^{p) ■= [ z{r)r'^dr > (m = 0, 1, 2). 
Jo 

Then we have 

Jlz{r){l-rfdr 2Z,{p) - Z2{p) 

z{r)dr Zq{p) 

Since < Z2{p) < Zi{p) holds, it is sufficient to show that Zi{p) / Zq{p) ^ as p — > 00. 
For m = 0, 1, we have the following equalities: 

Zm{p) = (y f e-P^^''^p{x)"'dx= [ [ eP(^(^)-^(^))(ixdy, 

Jo J J Dm 
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where 

Do := {{x, y) \ < X <l, <y <l}, Di := {(x, y) | < y < x < 1}. 

We define 

ujm'-= max {h{y) — b(x)) (m = 0, 1). 

{x,y)eDm 

Then, loq = bo and Assumption 15.11 is equivalent to ujq > uji. Defining 

G := y)eDo\ b{y) - b{x) > ^^^^ 

we obtain |G| > and 

Zi{p) _ Hp, ePWyy^i^y-^Uxdy ^ JJ^^ dxdy _ \D,\ 

as p tends to infinity. Hence we have proved (jS.lOp . □ 

As we have seen in the above proof, the exponential decay phenomenon of the principal 
eigenvalue occurs under the Assumption 15.11 

Corollary 5.3. Under A s sumption \ 5. 1\ we have 

lim -log = bo. (5.12) 

In other words, for arbitrary e £ (0, bo), there exists po > such that 

g-pCf-o+e) < ^^(p) < g-p(feo-e) (p>po). (5.13) 

Proof. It is sufficient to show (j5.12p , since ()5.13p is equivalent to (|5.12p . From Theorem 15.21 
Xi{p)Zo{p) ^ 1 as p ^ oo holds. We have (15.121) by the equalities: 

lim (- log -1- - bo] = lim (- logZo(p) - = lim - log (e-P^°Zo(p)) = 0, (5.14) 

where the last equality is shown as follows. From the inequality e^^^^' Zo{p) < \Do\, it follows 
that. 



lim sup - log e-P'^'Zoip) < lim ^' = 0. 

p^oo p \ ) p^oo p 

Moreover, for arbitrary e > 0, we define Dq := {{x, y) G Do; b{y) — b{x) >bo — e]. Then, from 
\Dq\ > and the inequality: 

e-P^''Zo{p)= [ [ eP(^(f)-^(^)-''o)dxdy > [ [ e-P'dxdy = e-P'\D'ol 
J J Do J Jd^ 



we obtain 



hminf-log (e-P'^Zoip)] > lim + = -e. 

p— >oo p \ ) p— »oo p 



Since e > is arbitrary, the last equality of ()5.14p follows. □ 
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For concrete examples of a{x), more precise asymptotic behaviom's of Xi{p) can be derived 
from Theorem 15.21 For example, if |iJi||i32| > then we get 

Ai(p) ^ IP HP | e~^^° asp^oo, 

since the following estimates hold: 

/ e-P^(^)dx ~ |5i|e-P^i asp ^00, if|Bi|>0, 
^0 

^ as p — > 00, if I-B2I > 0. 

Jo 

In case that |i?i||i32| = 0, the following lemma is useful. 

Lemma 5.4. Let L > and /i > 0. Suppose that g G C^{[0, L]) satisfies g{x) > for x G (0, L] 

and limxiox~^g{x) = 1. Then we have 

J e"^^(^)(ix ~ r ^- + 1^ asp ^00, 

where F stands for the Gamma function F(C) = s''~^e~^ds. 

Proof. In case that g{x) = x^, by the change of variables s = px^, we have 



1 r 1 



^, , dx = — — ; — r / si^ e ds = — — / e ds — > 1, 

F(i//i + i)yo F(i/;u + i)^p^yo r(i/^)yo 

as p tends to infinity. Hence, the assertion follows for g{x) = x^. For general g{x), we omit a 
proof but the same asymptotic behaviour can be shown. □ 

By using this lemma, for example, we obtain the following proposition. 
Proposition 5.5. For the case (15. 2p with a > 1, we have 



1 

;a— 1 



r (i + 1) ^ 

In particular, if a{x) = x ( i.e., —u"{x) +pxu'{x) = \u{x) {—I < x < I)), then 

[2 3 

)^i{p) ^ y — I p^ e~~P as p ^ CO. (5.15) 
Proof. From Lemma 15.41 we have 

l\-pK-)dx= ['e-^-"dx^T(- + l](^Y" asp ^00. 
Jo Jo V« / 

Using the change of variables y = I — x, we also have 
t eP''^'=Ux 
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fe^-''dx = e^'" f e~'^-'P3{y)dyr.e^'''T{2){l^-^p)-^ as p ^ 00, 
Jo Jo 
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where g{y) := a I "(Z" — (/ — y)") satisfies the condition of Lemma 15.41 with /i = 1. From 



Theorem 



we obtain the assertion. 



□ 



Similarly, we can calculate that, for (|5.3p : 



f V2sml 



p2 e 



-(1— cos l)p 



TT 



(if < Z < vr) 

(if / = vr) 

(if TT < Z < 2tt) 



as p 



oo, 



and for ([51 



as p 



oo, if ; > V2. 



We remark that similar precise asymptotic expansion has been studied in de Groen [3] by 
means of precise approximation of the eigenfunctions under a different assumption. Our result 
(j5.15p coincides to the result obtained there ((8.11) in [4J). 

At the end of this section, we remark that Theorem l5.2l can not be straightforwardly extended 
to multi-dimensional cases. To see this, we assume that Q and a{x) satisfy the conditions of 
Proposition 13.81 with Xj = and that each ai{x) satisfies Assumption 15.11 Then, since 

K^) = Z]r=i biixi) with 6'-(xi) = ai{xi), we have 



X 



dy 



4-"J|Af^(p) asp^oo, 



i=l 



whereas, from Proposition 13.81 we obtain 

n 

Xi{p) = Y.X?{p). 



In this case, the asymptotic behaviour of Xi{p) is determined by the slowest decaying X^^\p). 
6 Chemotaxis effect in biological colonies 

A simple biological interpretation of the exponential decay phenomenon of principal eigenvalue 
is discussed in this section. We consider the following simple PDE model for aggregation phe- 
nomena in biology. 

Let v{x,t) > be a population density at a; G of an organism (e.g. bacteria), which 
exhibits an aggregation behaviour due to an attractive chemical substance. This property is 
called chemotaxis and the attractive chemical substance is called chemotaxis substance. We 
denote by c{x, t) > the density of the chemotaxis substance. Then one of the standard PDE 
description of the chemotaxis effect is given by the following equation. 

vt = Av — pdiv{v'Vc) inr2x(0,oo), (6.1) 

where p > stands for a chemotaxis parameter. If p is larger, v exhibits stronger chemotaxis. 
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Keller and Segel ([10], [IT]) and other great number of literatures (for example Diaz and 
Nagai [6j, Stevens [16] etc.) have proposed (|6.ip nonlinearly coupled with the production of 
chemotaxis substance by v: 

act = Ac — JC + V in J7 X (0, oo), 

with a suitable boundary condition for c, where a and 7 are nonnegative constants. There are 
several results on finite time blow-up of these systems if the initial density v{x, 0) is large enough. 
The blowup of the solution in the Keller-Segel model has attracted many mathematicians' in- 
terests and its mathematical mechanism has been studied. Please see the above literature and 
references therein for more details. Especially, a good short review for mathematical results on 
the Keller-Segel model is found in § 8 of ^16j . 

On the other hand, in under the zero Dirichlet boundary condition for v with a = 0, 
they proved that v decays exponentially to zero as t ^ 00 if n = 1 or if v{x,0) is sufficiently 
small for n > 2. The zero Dirichlet boundary condition for v biologically means the extinction 
of the organism at the boundary. The result in [6j can be interpreted as that, the colony formed 
by the chemotaxis aggregation can not sustain eternally against the diffusion and the boundary 
extinction effect (t; = on 90), if the colony size is not large enough. 

Let us consider a simpler linear version of this chemotaxis aggregation vs boundary extinction 
problem. We suppose that the chemotaxis substance c is a priori given and does not change in 
time, i.e. c = c{x), and we consider the initial-boundary value problem: 

vt = A.V — pdiv{vVc{x)) inr2x(0,cxD) 
< v = onaJlx(0,oo) (6.2) 

, v{-,0) = vo{x) in n, 

where vo{x) > is a given initial density function with vq ^ 0. It is mathematically well-known 
that the boundary extinction is always stronger than the chemotaxis aggregation. It decays 
exponentially as 

v{x,t) Ce-^^^'P^*vi{x,p) (t^oo), 

where Xi{p) is nothing but Xi{p) in Section H] with b{x) = —c{x), and vi{x,p) is the principal 
eigenfunction of iC*^. If b{x) = —c{x) has the potential well (Definition 14. 6p of depth a; > 0, 
\i{p) exhibits the exponential decay phenomenon: \i{p) ~ 0{e~^P) for large p. The potential 
well condition for —c{x) means that there is a concentration peak of the chemotaxis substance 
of hight uj. 

The reciprocal of the principal eigenvalue 1/Ai(p) stands for the life span of the colony, 
since the half-life of exponential decaying quantity e~^* is given by (log2)/A. Theorem 14.71 is 
interpreted as that if c{x) has a concentration peak of hight uj > then the life span of the 
colony is exponentially long as 

-4t « Ce'^P as p^oo. 
Ai(p) 

In other words, the colony can not sustain eternally, but it sustains in enough long time 0{e^P) 
if the chemotaxis parameter is sufficiently large. 

In Figure [23l a typical profile of colony is shown. This is the eigenfunction vi of (13.30 
with same a(x) (= — Vc(a;)) as Example 12.21 and p = 40. Although the hight of this colony is 
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vl 




Figure 23: Eigenfunction vi of ()3.3p corre- 
sponding to ui of Figure [TOl 




Figure 24: Eigenfunction wi of (|3.3p corresponding to ui of Figure [2T] (left), and a scaled figure 
enlarged in the vertical axis (right). 



decreasing in time according to the linear law, but the decay is exponentially slow for large p 
and it substantially forms a very solid colony. Another profile of v corresponding to Example 12.31 
with p = 100 is shown in Figure [Ml which has two colonies. The hight of the right small colony 
is approximately 1/500 of the hight of the left one. This is interesting in comparison with 
Figure [2T] where the eigenfunction ui has two terraces with almost same hight. 
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